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Abstract 

In this talk I discuss the form factor approach used to compute correlation 
functions of integrable models in two dimensions. The Sinh-Gordon model 
is our basic example. Using Watson's and the recursive equations satisfied 
by matrix elements of local operators, I present the computation of the form 
factors of the elementary field <ft(x) and the stress-energy tensor T„ v (x) of the 
theory. 



1 Introduction 



Important progress in Quantum Field Theories has been made in recent years by 
studying two dimensional relativistic models. The reason essentially lays in the 
successful application of several non-perturbative approaches which lead to exact 
solutions of the quantum field dynamics. In this talk I am concerned with the 
determination of multi-point correlation functions of local operators for a massive 
field theory 

G n (n,r 2 ,...,r ft ) =<O il (r l )O i2 (r 2 )...O in (r n )> . (1.1) 

My basic example will the the Sinh-Gordon model which is discussed in a work 
made in collaboration with A. Fring and P. Simonetti [|J. 

For two-dimensional systems, there is usually a special limit where the compu- 
tation of the GnS is greatly simplified. This is the ultraviolet asymptotic regime of 
the correlation functions, defined by 

\ r i ~ r j\ < f • 

where the correlation length is given by £ = mT l . This regime can be equivalently 
reached by keeping finite values of the distances but taking the limit £ — > oo. Either 
ways, a scale invariant behaviour occurs in the model and the resulting 2-D infinite- 
dimensional conformal symmetry induces a system of linear differential equations 
satisfied by the correlators G n |§ . Their solution provides an explicit expression for 
the correlation functions of the theory in the ultraviolet regime 0, |3|, || . However, 
the behaviour of correlation functions in the cross-over region 

\n - r j\ - € , 

and in the infrared limit 
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is more complicated. Consider, for instance, the case of two-point correlation func- 
tions 

< Oi(r) 00) > = \ Hmr) • (1.2) 

With this parametrization, the exponent r)i is usually identified with twice the 
anomalous dimension of the field Oi whereas <&i(x) is a scaling function^] of the 
variable x = mr. Although conformal invariance severely restricts the possible 
values of r)i and gives quite powerful classification of the ultraviolet behaviours of 
two-dimensional models, the computation of the scaling functions $j(r) is, on the 
contrary, quite difficult. 

There exists, however, a large class of relativistic models where the determina- 
tion of the scaling functions can be explicitly worked out. This is the case of the 
massive integrable systems whose dynamics is strongly constrained by an infinite 
number of integrals of motion. In particular, the S'-matrix of these models presents 
factorization and elasticity properties and can be explicitly constructed [5-11]. In 
order to compute multi-point correlators for massive integrable models, we may 
exploit their spectral representation, i.e. their decomposition into an infinite sum 
over intermediate multi-particle contributions (fig. 1). For instance, the two-point 
function of a hermitian scalar operator Oi{x) in real Euclidean space can be written 
asQ 

(Oi{ X ) 0,(O)> = E / dfh U0 'tn n < °l<^ Wl' ...,&>**<&,..., Pn\O t (0)\0 > 

^ J n\(2Ti) n 

= E / I F n'(Pi ■ ■ -Pn) I 2 exp (-mr± cosh fl) , (1.4) 
go^ n!(27r)" V £i / 

* Scaling functions can be introduced as well for the n-point correlators. We concentrate our 

discussion on the two-point functions for clarity and simplicity. 

t/3 is the rapidity variable, related to the two-dimensional momentum by 

p° = m cosh/3 , p 1 — msinh/3 . (1.3) 
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where r denotes the radial distance r = Jxq + x\. 




Figure 1 

The functions 

F^(P 1 ,P 2 ,...,P n ) =<0\O i (0)\P 1 ,(h,...,Pn> i n (1-5) 

are the form factors of the operator 0, L (fig. 2). 




Figure 2 

Analytic properties of the form factors have been investigated by several authors 
and particularly important contributions can be found in [13-19]. As I will show in 
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the following, the form factor approach is a successful method in order to compute 
correlation functions. The expansion of the correlators in terms of the number of 
intermediate particles rather than in a perturbative series of the coupling constant 
presents, in fact, several advantages. First of all, the form factors take into account 
coupling constant dependence of correlation functions to all order in perturbation 
theory. Secondly, expressions like that one in eq. ( p..4|) are fast convergent series 
in the number N of intermediate particles even for moderate value of the distance 
and therefore, for practical applications, it is often sufficient to compute the form 
factors of the operator Oi with few external particles. Moreover, this spectral 
representation of two-point functions formally coincides with the grand partition 



function of a one-dimensional gas system, as noticed in [17, [LJ|. Hence, one can 



borrow successful techniques developped for one-dimensional system in order to 
extract non-perturbative parameters of the theory. For instance, as shown in [|T7], 



lQfl , the anomalous dimension of the operator O in the ultraviolet limit is equal to 



the pressure of the corresponding one-dimensional gas at one specific value of the 
fugacity. 

I will focalize on the form factor approach for the Sinh-Gordon model, described 
by the Lagrangian 

1 Tft ^ 

£ = ~(d^) 2 - -icoshg^x) . (1.6) 
^ 9 

The SG model belongs to the class of Toda Field Theory constructed on the simply- 
laced root systems, in this case A%. The reason why I will consider this model is 
its simplicity which allows us to understand the basic principles of the computation 
without being masked by algebraic complicancies. Let us discuss initially the basic 
features of this theory. 
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2 The Sinh-Gordon Theory 



The SG model is the simplest example of an affine Toda Field Theories |[20|| , pos- 
sessing a Zi symmetry — > —0. By an analytic continuation in g, i.e g — > ig>, it 
can formally be mapped to the Sine-Gordon model. 

There are numerous alternative viewpoints for the Sinh-Gordon model. First, 
it can be regarded either as a perturbation of the free massless conformal action by 
means of the relevant operator^ cosh g<f)(x). Alternatively, it can be considered as 
a perturbation of the conformal Liouville action 

1 



d x 



(2.1) 



by means of the relevant operator e~ 9 ^ or as a conformal affine Ai-Toda Theory 
pT[| in which the conformal symmetry is broken by setting the free field to zero. 

In a perturbative approach to the quantum field theory defined by the La- 
grangian ( |1.6|) , the only ultraviolet divergences which occur in any order in g come 
from tadpole graphs and can be removed by a normal ordering prescription with 
respect to an arbitrary mass scale M. All other Feynman graphs are convergent and 
give rise to finite wave function and mass renormalisation. The coupling constant 
g does not renormalise. 

An essential feature of the Sinh-Gordon theory is its integrability, which in the 



classical case can be established by means of the inverse scattering method [(24j] . 
In order to obtain the expressions of the (classical) conserved currents one can use 
Backhand transformation, as shown in jl], and obtains an infinite set of conservation 
laws 

d z T s+1 = ^6 S _! . (2.2) 



* Although the anomalous dimension of this operator (computed with respect to the free confor- 
mal point), is negative, A = — g 2 /8tt, the resulting theory is unitary. This is due to the existence 
of non a nonzero vacuum expectation values of some of the fields Oi in the theory. A detailed 
discussion of this point can be found in Q . 
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The corresponding charges Q s are given by 

Q* = <f [T s +i dz + 9 S _! dz] . (2.3) 

The integer-valued index s which labels the integrals of motion is the spin of the 
operators. Non-trivial conservation laws are obtained for odd values of s 

S =l,3,5,7,... (2.4) 

In analogy to the Sine-Gordon theory J23|, an infinite set of conserved charges Q s 



with spin s given in (2.4) also exists for the quantized version of the Sinh-Gordon 



theory. They are diagonalised by the asymptotic states with eigenvalues given by 

n 

Q s \P 1 ,...,P n >= XsE^'IA, •••,/?„> , (2.5) 



i=i 



where Xs is the normalization constant of the charge Q s . The existence of these 
higher integrals of motion precludes the possibility of production processes and 
hence guarantees that the n-particle scattering amplitudes are purely elastic and 
factorized into n(n — l)/2 two-particle S'-matrices. The exact expression for the 
Sinh-Gordon theory is given by |10| 



SM- ^-'Sl (2-6) 



tanh§03 + ^f ) 
where B is the following function of the coupling constant g 



b <«> = A- ■ (27) 



This formula has been checked against perturbation theory in ref. [[T(J and can also 
be obtained by analytic continuation of the S"-matrix of the first breather of the 
Sine-Gordon theory ||. For real values of g the S'-matrix has no poles in the physical 
sheet and hence there are no bound states, whereas two zeros are present at the 
crossing symmetric positions 

i— 

(3=1 2 (2.8) 



An interesting feature of the S'-matrix is its invariance under the map [|TTJ 



B^2-B (2.9) 
i.e. under the strong-weak coupling constant duality 

(2.10) 
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This duality is a property shared by the unperturbed conformal Liouville theory 
(|2.1|) p2]| and it is quite remarkable that it survives even when the conformal sym- 
metry is broken. 

3 Form Factors 

The form factors (FF) are matrix elements of local operators between the vacuum 
and n-particle in-state 

F°{p x ,p 2 ,...,p n ) =<0|O(0)|/3 1 ,/3 2 ,...,/3 n > in . (3.1) 

For local scalar operators 0(x), relativistic invariance implies that F n are functions 
of the difference of the rapidities. Except for the poles corresponding to the one- 
particle bound states in all sub-channels, we expect the form factors F n to be 
analytic inside the strip < Im/3jj < 2tt. 

The form factors of a hermitian local scalar operator 0(x) satisfy a set of func- 



tional equations, known as Watson's equations p2| , which for integrable systems 
assume a particularly simple form 

F n (fii, . . ■ , Pi, Pi+l, ■ ■ • , Pn) — F n (Px, . . . , Pi+l, Pi, . . . , Pn)S\Pi — Pi+l) , (3-2) 

n 

F n (p 1 + 2m,...,p n _ 1 ,p n ) = F n (p 2 ,...,p n ,p 1 ) = l[S(p i -P 1 )F n (p 1 ,...,P n ) . 

i=2 

The first equation states that as a result of the commutation of two particles in the 
asymptotic state we get a scattering process whereas the second equation fixes the 



discontinuity of the functions F n on the cuts fin = 2ni. In the case n = 2, eqs. 
reduce to 

W) = F 2 (-P)S 2 (/3) , 

F 2 (l7T-P) = F 2 (z7T + /3) . 

The general solution of Watson's equations for diagonal S-matrix systems can al- 
ways be brought into the form [T3| 



(3.3) 



< Prt 



K n ((3 1 ,...,(3 n )Y[F min ((3 l ,) , 

i<j 



(3.4) 



where F min ((3) has the properties that it satisfies ( |3.3j ), is analytic in < Im (3 < n, 
has no zeros in < Im (3 < it, and converges to a constant value for large values of 
f3. These requirements uniquely determine this function, up to a normalization. In 
the case of the SG model, F m - m {(5) is given by 



F min ((3,B) = J\f exp 



dx smh fef) sinh (f (1 - f )) sinh f f x p 
— ^5 sin — 



o x 



sinh x 



2tt 



We choose our normalization to be 



M 



exp 



rfx sinh(f ) S inh(f(l-f)) sinhf 



o x 



(3.5) 



(3.6) 



sinh x 

The analytic structure of F min (/3, B) can be easily read from its infinite product 
representation in terms of V functions 



f^(j3,b) = n 

k=0 



r(* + f + i)r(* + § + ! + |)r(* + i-f + | 



r(^ + | + H)r(A; + l-f + f)r(A; + i + f + f 

FminiP, B) has a simple zero at the threshold (3 = since S"(0) = — 
asymptotic behaviour is given by 

Jim F min (P,B) = 1. 

It satisfies the functional equation 

F min (in + P,B)F min (P,B) 



(3.7) 



1 and its 



sinh f3 



sinh p + sinh ^ 



(3.8) 



(3.9) 



S 



which is useful in order to find a convenient form for the recursive equations of the 
form factors. 

A useful expression for the numerical evaluation of F m j n (/3, B) is given by 



N-l 

F min (P,B) = A/"]} 

fc=0 



1 + 



1 + 



/3/27T 
K+2 4 



1 + 



( 0/27T 

U+i+f 



2 -te: 



/9/27r 



fe+i 



(3.10) 



x exp 



rf x sinh(f ) sinh(|(l-f)) sinhf 



(iV+l-iVe' 



•-' ) < 2Nx sin 2 I — 

io x sinh x \2tt / 

The rate of convergence of the integral may be improved substantially by increasing 
the value of N. 

The remaining factors K n in (|3.4j ) then satisfy Watson's equations with 5*2 = 1, 
which implies that they are completely symmetric, 27rz-periodic functions of the 
(3i. They must contain all the physical poles expected in the form factor under 
consideration and must satisfy a correct asymptotic behaviour for large value of 
Both requirements depend on the nature of the theory and on the operator O. 

Taking into account the one-particle pole in the three-particle channel at /3ij = 
in, the general form factors of a scalar hermitian operator in the Sinh-Gordon model 
can be parameterized as 



Fn(Pu ■ ■ ■ > Pn) — H n Q n (xx, 
where we have introduced the variables 



on 

Kj 



[Xi + X 



3, 



3 ft 



(3.11) 



(3.12) 



and the normalization constant H n . The functions Q n (%i, ■ ■ ■ ,%n) are symmetric 
polynomials^] in the variables £j. They can be expressed in terms of elementary 
symmetric polynomial a^ 1 \x±, . . . , x n ) which are generated by p5| 



W(x + Xi 

i=i 



J2 x 

k=0 



n-k (n)( 

°k \ X 1) X 2, • • • 5 X nJ 



(3.13) 



§The polynomial nature of the functions Q n is dictated by the locality of the theory |Tj 



Relativistic invariance demands that the total degree of Q n entering a form 
factor of spinless operators should be n(n — l)/2 (in order to match the total degree 
of the denominator in ( |3.11| )). The order of the degree of Q n in each variable Xi 
is fixed, on the other hand, by the nature and by the asymptotic behaviour of the 
operator O which is considered. We will come back to this point later, when we 
discuss specific operators. 

3.1 Pole Structure and Residue Equations for the Form 
Factors 

The pole structure of the form factors induces a set of recursive equations for the 
F n . For the Sinh-Gordon model there is only one kind of poles that arises from 
kinematical poles located at /3y = iir. They are related to the one-particle pole 
in a subchannel of three-particle states which, in turn, corresponds to a crossing 
process of the elastic S'-matrix. The corresponding residues are computed by the 
LSZ reduction [|l^, and give rise to a recursive equation between the n-particle 
and the (n + 2)-particle form factors (fig. 3) 

- i lim - (3)F n+2 + in, (3, ft, fo, ...,&)= ( 1 - f[ S(P - ft) J F n (P 1 , ...,&). 

/m* V i=i J 

(3.14) 
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This equation establishes a recursive structure between the (n + 2)- and n-particle 
form factors. 



4 Solution of the Recursive Equations 

According to the Z 2 symmetry of the SG model, we can label the operators by their 
parity. For operators which are Z 2 -odd the only possible non-zero form factors are 
those involving an odd number of particles. For Z 2 -even operators the only possible 
non-zero form factors are those involving an even number of particles. The vacuum 
expectation value of Z 2 -even operators can in principle be different from zero. 

The simplest representative of the odd sector is given by the (renormalised) field 
4>(x) itself. It creates a one-particle state from the vacuum. Our normalization is 
fixed to be 

Ff{p) =<O|0(O)|/3> in = i (4.1) 
For the even sector, an important operator is given by the energy-momentum tensor 

T^{x) = 2vr (: d^d v <j) - g^C{x) :) (4.2) 

where : : denotes the usual normal ordering prescription with respect an arbitrary 
mass scale M. Its trace TMx) = 0(x) is a spinless operator whose normalization 
is fixed in terms of its two-particle form factor 

F 2 @ ((3 12 = in) = out < fa | 6(0) | (3 2 > in = 2nm 2 , (4.3) 

where m is the physical mass. 

In the following we shall compute the form factors of the operators 4>{x) and 
Q(x). Employing the parameterization ( 3.11 ), the recursive equations ( 3.14Q take 



on the form 

(-) n Q n+2 (-x, x,xi,...,x n ) = xD n (x, xi,x 2 ,..., x n ) Q n (xi, x 2 ,..., x n ) (4.4) 
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where we have introduced the function 



D " = A ■ ( RlO\ 11 \( X + VXi)(x - U) 1 X i )\ - Yl \{X - UJXi)(x + LO 1 X i ) 

4sm^7ri5/ 1) \ i=1 i=1 ) 

(4.5) 

with u = exp(iTrB /2). The normalization constants for the form factors of odd and 
even operators are conveniently chosen to be 



F min (i7T, B) J 

'4Bin(7rS/2)Y 



n-l 

/ ! Sllll -/> / 7 1 \ 

H 2n = H 2 



where Hi and H 2 are the initial conditions, fixed by the nature of the operator. 
Using the generating function ( 3.13|) of the symmetric polynomials, the function 



D n can be expressed as 

D n = - . * f (-1Y sin ((k - l)^) x^a^a^. (4.7) 
2 sm(7rB/2) { ~ \ 2 J 

As function of B, D n is invariant under B — > —B. 

As shown in 0, the symmetric polynomials Q 2n +i entering the form factors of 
the elementary field 4>(x) can be factorized as 

Q2n+i(xi, ■ ■ ■ , x 2n+ i) = <J 2n +i "* P 2n +i(xi, ■ ■ ■ , x 2n+ i) n > , (4-8) 

whereas the analogous polynomials entering the form factors of the trace of the 
stress-energy tensor can be written as 

Q 2n {xi, . . . ,x 2n ) = a^ n) a^} 1 P 2n (x 1 , . . . ,x 2n ) n>l. (4.9) 

P n (xi, . . . ,x n ) is a symmetric polynomial of total degree n(n — 3)/2 and of degree 
n — 3 in each variable Xi. Using the following property of the elementary symmetric 
polynomials 

a k {-x,x,x 1 ,...,x n ) = a^\xi,x 2 , . . . ,x n ) - x 2 af] 2 {xi,x 2 , . . . ,x n ) , (4.10) 
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the recursive equations ( |4.4j ) can then be written in terms of the P n as 

( ) Pn+zi *^) "^i "^l; • • • ? •^n) D n (x, X±j X 2) ■ ■ ■ , 3^n) PnO^lj 2'2? • • • > ^<n) • 

(4.11) 

Using the recursive equations (|4.11|) and the transformation property of the elemen- 



tary symmetric polynomials (|4.10|) , the explicit expressions of the first polynomials 
P n (x u . . . , x n ) are given byfl 

P 3 (a;i, . . . ,x 3 ) = 1 
P 4 (xi, • • • ,x 4 ) = o 2 

P 5 (xi,..., x 5 ) = a 2 a 3 - cla 5 (4.12) 
P 6 (xi, ...,x 6 ) = o 3 (o 2 o 4 - cr 6 ) - c\(a 4 a b + a x a 2 a % ~ cr 3 cr 6 ) 
P 7 (xi, . . . , x 7 ) = o 2 o 3 o^ - cl(a 4 crl + aia 2 a 5 a 6 + ola 3 a 7 - c\a 2 a^a 7 ) + 
-02(0-10-20-40-7 + o 3 o 5 o 6 - c 2 oio 6 o 7 ) + c\c\a 7 

where c\ = 2cos(7rP/2) and c 2 = 1 — cf . Expression of the higher P n are eas- 
ily computed by an iterative use of eqs. Q4.4|) . For practical application the first 



representatives of P n are sufficient to compute with a high degree of accuracy the 
correlation functions of the fields. In fact, the n-particle term appearing in the 
correlation function of the fields ( |1.4| ) behaves as e _n(mr ) and for quite large values 
of mr the correlator is dominated by the lowest number of particle terms. This 
conclusion is also confirmed by an application of the c-theorem which I discuss at 
the end of the talk. It is interesting to notice that closed expressions for P n can be 
found for particular values of the coupling constant. 

4.0.1 The Self-Dual Point 



iThe upper index of the elementary symmetric polynomials entering P n is equal to n and we 
suppress it, in order to simplify the notation. 
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The self-dual point in the coupling constant manifold has the special value 

B (V8^) = 1 . (4.13) 



The two zeros of the ^-matrix merge together and the function D n (x, X\, x 2 , ■ ■ ■ , x n ) 
acquires the particularly simple form 



kir 



D n (x\x 1 ,x 2 ...,x n ) = ^(-l) fc+1 sin — 



X o 



kAn) 



£(-!)< COS ^"V^ 



\k=0 



U=0 



(4.14) 



In this case the general solution of the recursive equations ( |4.11 ) is given by [|TJ 

P n (xi,x 2 , . . . ,x n ) = detA(x 1 ,x 2 , . . . ,x n ) (4.15) 
where A is an (n — 3) x (n — 3) matrix whose entries are 



A ( \ _ ( n ) 2 

^ij\X\i %2i ■ ■ ■ 1 %n) — &2j— i+1 COS 



i.e. 



A 



( a 2 cr 6 • • • ^ 

a 3 cr 7 • • • 

1 0-4 ■■• 

0"i 05 • • • 



V 



/ 



This can be proved by exploiting the properties of determinants. 



(4.16) 



(4.17) 



4.0.2 The "Inverse Yang-Lee" Point 



A closed solution of the recursive equations (|4.11|) is also obtained for 

B(2^)=\. (4.18) 

The reason is that, for this particular value of the coupling constant the S-matrix 
of the Sinh-Gordon theory coincides with the inverse of the S'-matrix Syl{{3) of the 
Yang-Lee model |§ or, equivalently 



S(f3,--) = S YL (f3) . 



(4.19) 
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Since the recursive equations ( [4.11 ) are invariant under B — > —B (see sect. 4.2), a 
solution is provided by the same combination of symmetric polynomials found for 
the Yang-Lee model |L6L |j~8|, i.e. 



P n (x 1 ,x 2 , ...,x n ) = det B(xi, x 2 ,..., x n ) (4.20) 
with the following entries of the (n — 3) x {n — 3)-matrix B 

B i3 = * 3j -2i+i ■ (4.21) 

5 Form factors and c-theorem 

The Sinh-Gordon model can be regarded as deformation of the free massless theory 
with central charge c = 1. This fixed point governs the ultraviolet behaviour of the 
model whereas the infrared behaviour corresponds to a massive field theory with 
central charge c = 0. Going from the short- to large-distances, the variation of 



the central charge is dictated by the c-theorem of Zamolodchikov [26]. An integral 



version of this theorem has been derived by Cardy |Z7 and related to the spectral 
representation of the two-point function of the trace of the stress-energy tensor in 

mm 



i.e. 

POO 

Ac = / dfjLCi(fj,) , (5.1) 

J 



where ci(/x) is given by 

6 1 

G(p 2 ) = I d 2 xe~ ip± < 0|e(z)e(0)|0 > 



C M = ^— lmG (p 2 = -^) > ( 5 - 2 ) 



Inserting a complete set of in-state into (|5.2|) , we can express the function Ci(/z) in 
terms of the form factors Ff n 

12^2, 1 rdfr....^ , 



<*</.) = J^J I *£(A,..., A.) I (5.3) 



X 



5 m smb. j3i) 5(^2 m cosh Pi — ft) . 
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B 


9l 

47T 


Ac( 2 ) 


1 

500 

1 
100 

1 

10 

3 

10 

2 

5 

1 
2 

2 

3 

7 
10 

4 
5 

1 


2 
999 

2 
199 

2 

19 

6 
17 

1 

2 

2 
3 

1 

14 
13 

4 
3 

2 


0.9999995 
0.9999878 
0.9989538 
0.9931954 
0.9897087 
0.9863354 
0.9815944 
0.9808312 
0.9789824 
0.9774634 



For the Sinh-Gordon theory Ac = 1 and it is interesting to study the convergence 
of this series increasing the number of intermediate particles. For the two-particle 
contribution, we have the following expression 

The numerical results for different values of the coupling constant g 2 /47r are listed 
in the table below It is evident that the sum rule is saturated by the two-particle 
form factor also for large values of the coupling constant. Hence, the expansion 
in the number of intermediate particles results in a fast convergent series, as it is 
confirmed by the computation of the next terms involving the form factor with four 
and six particles. 
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6 Conclusions 



The computation of the Green functions is a central problem in a Quantum Field 
Theory. For integrable models, a promising approach to this question is given 
by the bootstrap principle applied to the computation of the matrix elements of 
local operators. It would be interesting to use this approach in order to derive 
differential equations satisfied by the quantum correlators and also to classify the 
operator content of a quantum integrable field theory. 



Acknowledgments 

I would like to thank A. Fring and P. Simonetti for our collaboration on this project 
and S. Elitzur and A. Schwimmer for useful discussions. 

References 

[1] A. Fring, G. Mussardo and P. Simonetti, Form Factors for Inte- 
grable Lagrangian Field Theories, the Sinh-Gordon Model, ISAS/EP/92-146, 
Imperial/TP/91-92/31, to appear on Nucl. Phys. B. 

[2] A. A. Belavin, A.M. Polyakov and A.B. Zamolodchikov, Nucl. Phys. B241 
(1984), 333. 

[3] Vl.S. Dotsenko and V.A. Fateev, Nucl. Phys. B240 [FS 12] (1984), 312; Nucl. 
Phys. B251 [FS 13] (1985), 691; Phys. Lett. B 154 (1985), 291. 

[4] C. Itzykson, H. Saleur and J.B. Zuber, Conformal Invariance and Applications 
to Statistical Mechanics, (World Scientific, Singapore 1988). 

[5] A.B. Zamolodchikov, ALB. Zamolodchikov, Ann.Phys. 120 (1979), 253. 

17 



[6] A.B. Zamolodchikov, in Advanced Studies in Pure Mathematics 19 (1989), 641; 
Int. J. Mod. Phys. A3 (1988), 743. 

[7] R. Kdberle and J.A. Swieca, Phys. Lett. 86B (1979), 209; A.B. Zamolodchikov, 
Int. J. Mod. Phys. A3 (1988), 743; V. A. Fateev and A.B. Zamolodchikov, Int. 
J. Mod. Phys. A5 (1990), 1025. 

[8] J.L. Cardy, G. Mussardo, Phys. Lett. B225 (1989) , 275. 

[9] G. Mussardo, Phys. Rep. 218 (1992), 215. 

[10] A.E. Arinshtein, V.A. Fateev and A.B. Zamolodchikov, Phys. Lett. 87B (1979), 
389. 

[11] P. Christe and G. Mussardo, Nucl.Phys. 5B330 (1990), 465; P. Christe and G. 
Mussardo Int. J. Mod. Phys. A5 (1990), 1025; H. W. Braden, E. Corrigan, P. 
E. Dorey, R. Sasaki, Nucl. Phys. B338 (1990), 689; H. W. Braden, E. Corrigan, 
P. E. Dorey, R. Sasaki, Nucl. Phys. B356 (1991), 469. 

[12] K.M. Watson, Phys. Rev. 95 (1954), 228. 

[13] B. Berg, M. Karowski, P. Weisz, Phys. Rev. D19 (1979), 2477; M. Karowski, 
P. Weisz, Nucl. Phys. B139 (1978), 445; M. Karowski, Phys. Rep. 49 (1979), 
229; 

[14] F. A. Smirnov, in Introduction to Quantum Group and Integrable Massive Mod- 
els of Quantum Field Theory, Nankai Lectures on Mathematical Physics, World 
Scientific 1990. 

[15] F.A. Smirnov, J. Phys. A17 (1984), L873; F.A. Smirnov, J. Phys. A19 (1984), 
L575; A.N. Kirillov and F.A. Smirnov, Phys. Lett. B198 (1987), 506; A.N. 
Kirillov and F.A. Smirnov, Int. J. Mod. Phys. A3 (1988), 731. 



18 



[16] F.A. Smirnov, Nucl. Phys. B337 (1989), 156; Int. J. Mod. Phys. A4 (1989), 
4213. 

[17] V.P. Yurov and Al. B. Zamolodchikov, Int. J. Mod. Phys. A6 (1991), 3419. 

[18] ALB. Zamolodchikov, Nucl. Phys. B348 (1991), 619. 

[19] J.L. Cardy and G. Mussardo, Nucl. Phys. B340 (1990), 387. 

[20] A.V. Mikhailov, M.A. Olshanetsky and A.M. Perelomov, Comm. Math. Phys. 
79 (1981), 473. 

[21] O. Babelon and L. Bonora, Phys. Lett. B244 (1990), 220. 
[22] P. Mansfield, Nucl. Phys. B222 (1983), 419. 

[23] R. Sasaki and I. Yamanaka, in Advanced Studies in Pure Mathematics 16 
(1988), 271. 

[24] L.D. Faddev and L.A. Takhtajan, Hamiltonian Method in the Theory of Soli- 
tons, (Springer, N.Y., 1987). 

[25] I.G. MacDonald, Symmetric Functions and Hall Polynomials (Clarendon Press, 
Oxford, 1979). 

[26] A.B. Zamolodchikov, JEPT Lett. 43 (1986), 730. 

[27] J.L. Cardy, Phys. Rev. Lett. 60 (1988), 2709. 

[28] A. Cappelli, D. Friedan and J.L. Latorre, Nucl. Phys. B352 (1991), 616. 
[29] D.Z. Freedman, J.L Latorre and X. Vilasis, Mod. Phys. Lett. A6 (1991), 531. 

19 



